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P1552 

SYSTEM AND METHOD FOR AUTOMATIC DETERMINATION OF 
A REGION OF INTEREST WITHIN AN IMAGE 

TECHNICAL FIELD 

[0001] The present invention relates to automated determination of a Region Of Interest 
(ROI) within an Image, and more particularly, to an automated determination of an ROI of an 
ultra sound image of a carotid artery. 

BACKGROUND ART 

[0002] According to the U. S. Center for Disease Control (CDC), heart disease is the leading 
cause of death in the United States, and atherosclerosis (a progressive narrowing and hardening 
of the arteries over time) is the leading factor in heart disease. Although initial thickening of 
arterial walls develops over a very long period of time, the complications from atherosclerosis 
generally occur suddenly and often prove fatal. Further, because atherosclerosis is asymptomatic, 
the first symptom perceived by a patient may be a life-threatening event. 
[0003] Several methods are available for measuring the progress of atherosclerosis. These 
include Electron Beam Computed Tomography (EBCT), Magnetic Resonance Imaging (MRI), 
Flow Mediated Dilation (FMD), Quantitative Intima Media Thickness (QIMT), and the like. 
Although MRI may be use for clinical measurements, it remains a costly procedure. FMD 
appears promising, but outcome studies are not available to validate its predictive ability. 
[0004] Carotid Intima Media Thickness (IMT) has been shown to be a useful surrogate 
endpoint. A surrogate endpoint is an indirect measurement which correlates with an actual (i.e., 
clinical) outcome. A clinical endpoint, for example, could be death, loss of vision, or other major 
symptomatic event. A corresponding surrogate endpoint could be a measurement which 
indicates the development of a condition which is correlated with a clinical endpoint event, but 
which surrogate endpoint may not produce easily observable symptoms. In particular, IMT 
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measurements of the far wall of the carotid artery, about 1 cm below the bulb, have shown high 
correlation with cardiovascular clinical endpoints. 

[0005] Carotid EMT provides a suitable surrogate endpoint for all stages of atherosclerosis 
and is well suited for tracking the progress of atherosclerosis. Carotid IMT has predictive value 
with respect to cerebrovascular and cardiovascular clinical events in both men and women. It 
has further been found that providing an ultrasound image of the carotid artery to a patient has 
a substantial effect of motivating changes to behavior, and the provision of subsequent images 
showing progress helps patients remain motivated to sustain a healthy life style. 
[0006] Carotid IMT may be measured by either analog caliper IMT or digital Quantitative 
IMT (QIMT). Both methods rely on the same image and attempt to measure the same quantity. 
Caliper IMT is a manual measurement of the IMT visible in an ultrasound image using calipers, 
and is similar to printing a photo, and then measuring a dimension of a characteristic of the 
photograph by laying a ruler on the photograph. QIMT is a computer analysis of IMT performed 
on a digital version of the ultrasound image. Results have shown that caliper IMT has four times 
the standard deviation of QIMT when performed on the same images. Because of the accuracy 
of QIMT measurements, QIMT is ideal for both absolute measurements of the IMT, and for 
tracking small changes in the carotid artery, and thus the progression of atherosclerosis. The 
ability to track small changes allows the success or failure of a given course of treatment to be 
analyzed and adjusted if necessary, before a critical state is reached. 

[0007] Current practice in the field of QIMT is to record carotid artery ultrasound images on 
SVHS tapes, and later manually review the images to determine a Region Of Interest (ROI) 
within each frame for subsequent computer analysis using edge detection algorithms. Such 
manual review requires highly skilled technicians and is time intensive. Because of the present 
requirement for manual review, IMT presently is only available from specialized laboratories 
employing highly trained technicians. There is thus a need for an objective system and method 
to automate the determination of the ROI of carotid ultrasound images, which method does not 
introduce operator bias. 

DISCLOSURE OF THE INVENTION 
[0008] The present invention addresses the above and other needs by providing a system and 
method for automated image analysis, which system and method selects a Region Of Interest 
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(ROI) within the image (or frame) suitable for subsequent Quantitative Intima Media Thickness 
(QIMT) analysis. Presently, a time-consuming manual analysis by a highly trained technician 
is required to select the ROI within the frame for subsequent analysis. The automated image 
analysis of the present invention provides an efficient, timely, and accurate identification of the 
ROI within the frame. 

[0009] In accordance with one aspect of the invention, there is provided a method for 
defining an ROI for subsequent processing, comprising specifying a physical characteristic of a 
morphological structure, providing at least one image of a portion of the morphological structure, 
said image containing the physical characteristic, and automatically processing the image to 
identify an ROI within the image, wherein the ROI contains the physical characteristic, and 
wherein the ROI is smaller than the original image, and the ROI is preferably from approximately 
0.25 centimeters to approximately 1.0 centimeters wide by approximately 0.25 centimeters to 
approximately 1.0 centimeters high, and the ROI is more preferably approximately 0.5 
centimeters wide by approximately 0.5 centimeters high. 

[0010] The physical characteristic is preferably an Intima Media Thickness (IMT) of an 
artery, more preferably the IMT of a portion of a carotid artery, and most preferably the IMT of 
a far wall of the Bulb-Common Carotid Artery interface ranging between approximately LO cm 
and approximately 1 .5 cm below the base of the bulb, or most preferably the IMT of the far wall 
of a portion of the carotid artery ranging between approximately 1 .75 cm and approximately 2.25 
cm below the bifurcation. The image is preferably processed digitally, and the image may be 
provided as a digital image, or the image may be provided as an analog image, and processed to 
generate a digital image, and more preferably, the image is provided as a SVHS image and is 
processed to generate a digital image. 

[0011] In one embodiment, the image is provided as a rectangular frame comprising a 
rectangular array of pixels, which frame is vertically sliced into a multiplicity of slices, which 
slices are preferably two pixels wide. The slices are processed to identify the presence of stripes. 
A stripe comprises three bands, a dark band with light bands both above and below the dark 
band, which light bands correspond to vessel walls, and which dark band corresponds to a vessel 
lumen. Each slice is processed by image processing software which identifies the stripes with 
confidence scores above a threshold, preferably a fifty percent confidence threshold. 
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[001 2] Pairs of slices may be compared to detect the bifurcation of the carotid artery, wherein 
the pairs of sUces preferably comprise adjacent slices. If for any reason, adjacent slices cannot 
be singly compared, substantially adjacent slices (e.g., less than 3 slices apart) maybe compared, 
or a combination of substantially adjacent slices maybe compared (i.e., an averaging of features). 
If a first slice (i.e., a slice above, or after another slice) of a pair of slices has one more stripe than 
a second slice (i.e., a slice below or before another slice) and if top most and bottom most edges 
of two adjacent first stripes in the first slice (i.e., the slice with a greater number of stripes) 
correspond respectively to the top edge and bottom edge of a second stripe in the second slice, 
a bifiircation is identified. Preferably, the vertical positions of the corresponding edges will be 
within approximately five percent of each other. If the bifiircation is foimd, the ROI is 
horizontally centered approximately 2.0 centimeters before the bifiircation. 
[0013] In some cases, the bifiircation may not be successfiiUy identified. In these cases, the 
slices may be processed to attempt to identify the bulb of the carotid artery. The average number 
of stripes within the fi*ame is computed. If a valid number of stripes is computed, an attempt is 
made to use the resulting value to determine the position of the bulb of the carotid artery. The 
portion of the carotid artery below the bulb is referred to as the Common Carotid Artery (CCA). 
If the number of stripes is three or more, the stripe that is next to the bottom is tagged as a selected 
vessel. If the number of stripes is two, the bottom stripe is tagged as the selected vessel. If the 
number of stripes is one, the one stripe is tagged as the selected vessel. 

[0014] A moving average of the vessel diameter (i.e., lumen diameter) is computed, preferably 
from ten consecutive slices. The center of the moving average of slices is tagged as the bulb 
position if the lumen diameter increases sufficiently fi"om one average to the next average, 
preferably at least ten percent lumen diameter increase. If the bulb is found, the ROI is 
horizontally centered approximately 1.25 centimeters before the bulb. 

[0015] In some cases, neither the bifurcation nor the bulb are detectable. In these cases, the 
automatic processing may select a defined region approximately twenty percent higher and 
approximately twenty percent wider than the ROI from a previous fi-ame. A pattern match may 
then be automatically performed within the defined region to determine a new ROI which matches 
a previous ROI. If all of the automatic processing fails to identify the ROI, the ROI may be a 
manually generated by a technician. 
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[0016] It is a further feature of the invention to provide a method for determining a Region 
Of Interest (ROI) for subsequent Intima Media Thickness (IMT) analysis, the method comprising 
providing at least one digital frame of a portion of a carotid artery, wherein the frame includes at 
least one of a set consisting of a bulb of the carotid artery and a bifrircation of the carotid artery, 
and automatically processing the frame to identify an ROI in the frame. The automatically 
processing comprises at least one of the set consisting of determining the ROI by detecting a 
bifurcation of the carotid artery, determining the ROI by detecting a bulb of the carotid artery, and 
determining the ROI by pattern-matching a portion of the frame to a previous ROI. 
[001 7] Determining the ROI by detecting a bifurcation of the carotid artery comprises slicing 
the frame into vertical slices, automatically comparing at least one of the vertical slices to at least 
one other of the vertical slices to find the bifurcation of the carotid artery, determining an ROI 
horizontal extent ranging from approximately 1.75 cm before the bifurcation to approximately 
2.25 cm before the bifurcation, identifying a highest point of a far wall segment of the carotid 
artery within the ROI horizontal extent, and determining an ROI vertical extent ranging from 
approximately 0.3 cm above the highest point to approximately 0.2 cm below the highest point. 
[0018] Determining the ROI by detecting a bulb of the carotid artery comprises slicing the 
frame into the vertical slices, automatically comparing at least one of the vertical slices to at least 
one other of the vertical slices to find the bulb of the carotid artery, determining the ROI 
horizontal extent to range from approximately 1 .00 cm before the bulb to approximately 1 .50 cm 
before the bulb, identifying the highest point of a far wall segment of the carotid artery within the 
ROI horizontal extent, and determining the ROI vertical extent to range from approximately 
0.3 cm above the highest point to approximately 0.2 cm below the highest point. 
[0019] Determining the ROI by pattem-matching a portion of the frame to a previous ROI 
having a previous horizontal extent and a previous vertical extent comprises determining an 
expanded ROI horizontal extent to the previous horizontal extent expanded horizontally by a 
factor of approximately 1 .2, determining an expanded ROI vertical extent to the previous vertical 
extent expanded vertically by a factor of approximately 1.2, performing an automatic pattern 
match comparing the previous ROI to the region defined by the expanded ROI horizontal extent 
and the expanded ROI vertical extent, and determining the ROI based on the highest correlation 
of the automatic pattem match. 
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[0020] It is an additional feature of the invention to provide a system for processing ultrasound 
images of the carotid artery, comprising an ultrasound sensor, an ultrasound apparatus adapted to 
generate frames, a sensor cable connecting the ultrasoimd sensor to the ultrasound apparatus, a 
computer, and a path adapted to carry the ultrasound data from the ultrasound apparatus to the 
computer, wherein the path is selected from a group consisting of dedicated cables, phone lines, 
a Local Area Network (LAN), the Internet, and a wireless path. 

[0021] The ultrasound apparatus is adapted to receive measurements from the ultrasound 
sensor and to process the measurements to generate ultrasound data. The computer is adapted to 
pre-process the ultrasound data to generate Regions Of Interest (ROI) within the data for 
subsequent processing, which ROIs are subsets of the ultrasound data. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0022] The above and other aspects, features and advantages of the present invention will be 
more apparent from the following more particular description thereof, presented in conjunction 
with the following drawings wherein: 

[0023] FIG. 1 shows a system suitable for acquiring and processing ultrasound image data 
according to the present invention; 

[0024] FIG. 2 depicts positioning of an ultrasoimd sensor relative to a carotid artery; 
[0025] FIG. 3 depicts a flow chart of automated Region Of Interest (ROI) determination; 
[0026] FIG. 4A shows an ultrasound image frame; 

[0027] FIG. 4B shows an altemative representation of an ultrasoimd image; 

[0028] FIG. 4C shows how the image of FIG. 4A may be sliced into slices for processing 

according to the method of the present invention; 

[0029] FIG. 5 shows an example of slices which may be compared to determine a position of 
a bifurcation of the carotid artery; 

[0030] FIG. 6 depicts the placement of the ROI relative to the bifurcation of the carotid artery; 
[0031] FIG. 7 depicts a single slice suitable for processing to determine a number of stripes 
(or vessels) in the slice; 

[0032] FIG. 8 depicts the placement of the ROI relative to a bulb of the carotid artery 
[0033] FIG. 9 depicts placement of an expanded region for a pattern match around a 
previously determined ROI; 
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[0034] FIG. 10 shows a flow chart of a method for analyzing a frame to find the bifurcation 
and/or the number of vessels within the frame; 

[0035] FIG. 1 1 shows a flow chart of a method for analyzing a frame to find an average 
number of vessels in the frame; and 

[0036] FIG. 12 shows a flow chart of a method for technician-generated ROI. 

[0037] Corresponding reference characters indicate corresponding components throughout the 

several views of the drawings. 

MODES FOR CARRYING OUT THE INVENTION 
[0038] The following description is of the best mode presently contemplated for carrying out 
the invention. This description is not to be taken in a limiting sense, but is made merely for the 
purpose of describing the general principles of the invention. The scope of the invention should 
be determined with reference to the claims. 

[0039] The present invention provides a system and method for automated image quality 
analysis which selects a Region Of Interest (ROI) within the frame for subsequent analysis. In one 
embodiment, the subsequent analysis comprises computing a Quantitative Mima Media 
Thickness (QIMT) of a carotid artery. The QIMT is a usefiil parameter in determining the 
cardiovascular health of a patient, and particularly the presence of heart disease, and more 
particularly the presence of atherosclerosis. QIMT is capable of providing very accurate 
measurements, and as a result QIMT is ideal for both absolute measurements of the IMT, and for 
tracking small changes in the carotid artery, and thus the progression of atherosclerosis. The 
ability to track small changes allows the success or failure of a given course of treatment to be 
analyzed and adjusted if necessary, before a critical state is reached. 

[0040] QIMT may be performed using, for example, the ARTIS computerized edge contour 
measurement technique developed by Prevention Concepts Inc., in West Los Angeles, Califomia. 
Although the automated image quaUty analysis method of the present invention is described herein 
in the context of pre-processing images for QIMT of the carotid artery, those skilled in the art will 
recognize that other methods of processing images will similarly benefit from the present 
invention, and in particular, other methods of processing ultrasound images will similarly benefit 
from the present invention, and even more particularly, other methods of processing ultrasound 
images of any artery or of the carotid artery in particular will similarly benefit from the present 
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invention, and the application of the present invention to pre-processing any image is intended to 
come within the scope of the present invention. 

[0041] Elements of a typical system according to the present invention is shown in FIG. 1 . 
The elements include an ultrasound sensor 12 shown residing proximal to the neck of a patient 
10. The ultrasound sensor 12 is connected by a sensor cable 14 to an ultrasound apparatus 16. 
The ultrasoimd apparatus 16 may provide power to the ultrasound sensor 12 and receive 
measurements from the ultrasound sensor 12. The ultrasound apparatus 16 communicates with 
a computer 20 over a path 18. The path 1 8 may be any of a variety of paths including dedicated 
cables, phone lines, a Local Area Network (LAN), the Intemet, a wireless path, etc., and a system 
according to the present invention using any path is intended to come within the scope of the 
present invention. 

[0042] The computer 20 processes the signal received from the ultrasound apparatus 16 and 
may provide results on a screen 26, or over an output cable 22 to an output device 24. When the 
computer 20 is remote from the ultrasound apparatus 16, the results may be transmitted back to 
the ultrasound apparatus over the path 18, or over a different path, and be displayed on a second 
screen 28 

[0043] A more detailed view of the patient 10 and ultrasound sensor 12 is shown in FIG. 2. 
The ultrasound sensor 12 is positioned to image a far wall 32 of the carotid artery 30, about one 
cm before a bulb 34, or about two centimeters before a bifiircation 35. The blood flow within the 
carotid artery 30 is from the right to the left in FIG. 2, and a point within the carotid artery 30 to 
the right of the bulb 34 or to the right of the bifiircation 35 is hereafter referred to as being before 
or below the bulb 34 or the bifiircation 35. Various studies have shown that measurements the 
IMT of the far wall 32 correlate well with clinical end points. The far wall 32 IMT measurements 
have proven to be reproducible, and show a remarkable correlation with cardiovascular events as 
far as ten years in the fiiture. 

[0044] A method for automated Region Of Interest (ROI) determination according to the 
present invention is shown in FIG. 3. Ultrasound image data of the carotid artery 30 is collected, 
for example as depicted in FIG. 2 above, in a format of image frames. The image frames may be 
provided as digital images, or the image frames may be provided as analog images and processed 
to generate digital images, and preferably, the image frames may be provided as SVHS images 
and processed to generate digital images. The image frames are provided to select frame for 
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analysis at 42 where a frame is selected for subsequent processing. The frame is analyzed for 
bifrircation and number of vessels at 44. If a bifurcation is found at 46, the ROI is placed based 
on the location of the bifurcation at 48. If a bifurcation is not found, but at least one vessel is 
found 50, the ROI is placed based on the vessel geometry at 52. If neither bifurcation nor at least 
one vessel is found, an attempt is made to place the ROI based on a pattem match at 54. If a 
match was foimd at 56, a new frame is selected for processing. If a match was not found, and if 
only a small number of frames have failed to support ROI placement, the ROI is placed as in the 
previous frame at 58. If more than a small number of frames have failed to support ROI 
placement, the ROI is generated manually by a technician at 62. 

[0045] The image frames are generally provided as rectangular first frame 86a as shown in 
FIG. 4A. Each frame 86a may include a multiplicity of exemplar bands 90a - 90g, and a 
bifurcation 91 of the carotid artery (FIG. 2). The bands 90a -90g are for purposes of illustration 
only. The bands in an actual frame may be dark or light, and may be narrow or wide. An example 
of a second frame 86b displayed to show vessels as broad gray bands is shown in FIG. 4B. The 
frames 86a and 86b preferably comprise a rectangular array of pixels. 

[0046] The frames 86a, 86b may be vertically sliced into a multiplicity of slices 88 as shown 
in FIG, 4C, which slices are preferably two pixels wide. A first slice 92a and a second slice 92b 
are labeled, which slices 92a, 92b straddle the bifurcation 91. The slices 92a, 92b are processed 
to identify the presence of stripes. A stripe comprises three bands, a dark band with Ught bands 
both above and below the dark band, which light bands correspond to vessel walls, and which 
dark band corresponds to a vessel lumen. Each slice is processed by image processing software 
which identifies the stripes and computes confidence scores, stripes used in subsequent processing 
preferably exceed a fifty percent confidence threshold. Commonly available edge detection 
software maybe used for processing the images, preferably the Active-Mil 7.0 toolkit developed 
by Matrox Electronic System Ltd. in Dorval, Quebec, Canada, or Halcon 6.1 developed by 
MVTec Software GmbH in Muchen, Germany. 

[0047] The slices maybe compared to detect the bifurcation of the carotid artery, preferably 
adjacent slices are compared. If for any reason, adjacent slices can not be singly compared, 
substantially adjacent slices (i.e., less than 3 slices away) maybe compared, or a combination of 
substantially adjacent slices may be compared (i.e., an averaging of features). For example, the 
first slice 92a and second slice 92b maybe compared as shown in FIG. 5, wherein Stripes SI, S2, 
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and S3 are identified in slice 92a, and stripes S4 and S5 are identified in slice 92b. If the number 
of stripes in slice 92a is one more than the number of stripes in slice 92b, and if the top edge stripe 
S2 corresponds to the top edge of stripe S5, and the bottom edge of stripe S3 corresponds to the 
bottom edge of stripe S5, a bifiircation 91 is identified. Preferably, the vertical positions of the 
corresponding edges will be within five percent of each other. 

[0048] If the bifiircation 91 is identified, the ROI is horizontally centered approximately 2.0 
centimeters before the bifiircation, shown by dimension Dl in FIG. 6. The ROI defines a top edge 
104, a bottom edge 102, a right edge 106, and a left edge 100. The ROI is preferably 
approximately 0.5 centimeters wide (W), and preferably approximately 0.5 centimeters high (H), 
and is vertically positioned with a top edge 104 approximately 0.3 centimeters above highest 
point of a bottom vessel wall (i.e., the far wall) within the horizontal range of the ROI, of the 
vessel containing the bifiircation. 

[0049] In some cases, the bifiircation may not be successfiiUy identified. In these cases, the 
slices may be processed to attempt to identify the bulb of the carotid artery. The number of stripes 
(or vessels) in each slice, for example in a third slice 92c as shown in FIG. 7, is computed. An 
average number of stripes in the fi-ame is computed by averaging the number of stripes in all of 
the slices within the fi-ame, preferably the average number of stripes in computed using stripes 
with high confidence, and more preferably the average number of stripes is computed using stripes 
having at least a fifty percent confidence score. 

[0050] A rounded average number of stripes is computed by roimding the average to the 
nearest whole number. The mean and variance of the number of stripes is computed relative to 
the rounded average number of stripes. The roimded average number of stripes is tagged as valid 
if the mean is below a mean threshold, and the variance is below a var threshold. The 
mean_threshold is preferably approximately 0.25, and the var_threshold is preferably 
approximately 1 .0. 

[0051] The minimum confidence level of the stripes within each slice is compared to a 
min_con_threshold. Slices having one or more stripes with a confidence level below the 
min_con_threshold are tagged as outliers, and not used in subsequent computations. A second 
variance is computed based on the remaining slices, and the second variance is compared to a 
second_var_threshold. An average confidence level is computed for the remaining stripes and 
compared to an avg_con_threshold. If either the second variance is greater than the 
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second_var_threshold, or the average confidence is less than the avg_con_threshold, the average 
number of stripes is tagged as invalid. The second_var_threshold is preferably approximately 1 ,0, 
and the avg_con_threshold is preferably approximately 75 percent. 

[0052] If a valid number of stripes is computed, average_stripes is set to the rounded average 
number of stripes, and an attempt is made to use the resulting value to determine the position of 
the bulb 34 of the carotid artery 30 (see FIG. 2). The portion of the carotid artery below the bulb 
is referred to as the Common Carotid Artery (CCA). If the average_stripes is three or more, the 
stripe that is next to the bottom is tagged as a selected stripe. If the average_stripes is two, the 
bottom stripe is tagged as the selected stripe. If the average_stripes is one, the one stripe in each 
slice with a single stripe is tagged as the selected stripe. 

[0053] In each of the three cases described above, one or more of the slices may have a 
number of stripes different from the average_stripes. In this case, initially, the stripes are selected 
from the slices having a number of stripes equal to the average_stripes. Stripe position is set to 
the average of the positions of the tagged stripes. The stripe in the remaining slices (i.e., slices 
with more than one stripe) which overlaps the stripe position is tagged as the selected stripe. If 
no stripe overlaps the stripe position, the slice is not used for subsequent processing. 
[0054] A moving average of the stripe diameters is computed using several consecutive 
selected stripes, the nxmiber of consecutive selected stripes is preferably between five and fifteen, 
and. the niunber of consecutive selected stripes is more preferably ten. Consecutive averages are 
compared. If the average increases by at least ten percent, the bulb position 112 is set to the 
position of the center slice from the average, as shown in FIG. 8. The RIO is horizontally centered 
approximately 1.25 centimeters before the bulb (shown as D2). The ROI is vertically centered 
0.05 centimeters above the highest point of the far wall of the carotid artery (i.e., the bottom of 
the selected stripe of the frame) within the horizontal range of the ROI. The ROI width (W) is 
preferably 0.5 centimeters and the ROI height (H) is preferably 0.5 centimeters. 
[0055] If neither the bifiircation nor the bulb are found, a pattem match may be attempted 
within a failed frame 1 14 to position the ROI as shown in FIG, 9. An expanded ROI horizontal 
extent is set to the previous horizontal extent expanded horizontally by a factor of approximately 
1 .2. An expanded ROI vertical extent is set to the previous vertical extent expanded vertically by 
a factor of approximately 1 .2. The ROI center is set to the center of the previous ROI. A pattem 
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match is then performed in the expanded region 1 1 6 to find an area that matches the previous ROI 
118. Preferably, pattern-matching is performed using the Matrox Toolkit, or the MVTec toolkit. 
[0056] If neither the bifurcation, nor the bulb, nor a pattern match are found, and if only a 
small number of consecutive failures have occurred, the ROI is set to the previous ROI, otherwise, 
the ROI is manually determined by a technician. 

[0057] A method for analyzing a frame for bifurcation 44 is shown in FIG. 10. A first vertical 
slice is selected for processing at 66, and first stripes in the first slice are identified at 68. A 
second vertical slice is selected for processing at 70, and second stripes within the second slice 
are identified at 72. The first and second stripes are compared to attempt to determine if the 
bifurcation occurs substantially between the first and the second slice at 74 (i.e., the bifurcation 
is in the region substantially overlapping the first and second slices). If the bifurcation is foimd 
at 76, the bifurcation location is saved at 78. Otherwise, if all slices have been compared at 82, 
the average number of slices (i.e., vessels) is computed and saved at 84. Otherwise, the 
processing starting with select a first vertical slice 66 is repeated for new first and/or second slices. 
[0058] A method for determining an average number of vessels in a frame 84 (see FIG. 10) 
is shown in FIG, 1 1 . An average number of stripes is computed at 86. A roimded average number 
of stripes is computed from the average number of stripes at 88. The mean and variance of the 
number of stripes is computed at 90. If the mean is not below a mean threshold M (preferably 
0.25) and the variance is not below a first variance threshold VI (preferably 1 .0) at 92, the average 
is tagged as not valid at 94. Otherwise, slices with one or more low confidence stripes are 
discarded at 96. A second average number of stripes is computed from remaining slices at 98. 
A second variance and average confidence are computed at 100. If the average confidence is 
below a threshold (preferably seventy five percent) and the second variance is above a second 
variance threshold (preferably 1.0) at 102, the average is tagged as not valid at 104. Otherwise, 
the number of vessels is saved for subsequent processing at 106. 

[0059] A method for technician-generated ROI 62 (see FIG. 3) is shown in FIG. 12. The ROI 
is determined by the technician at 129. A new frame is selected for analysis at 122. A pattern 
match is attempted at 124, which match is similar to the match described in FIG. 9 above. If a 
match is found at 1 26, the technician may reinitialize automatic ROI determination. If automatic 
ROI determination is reinitialized at 136, automode is returned to, otherwise another new frame 
is selected at 122 and the processing described above is repeated. 
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[0060] If a match was not found at 1 26, and more than a small number of frames have failed 
to result in a matching pattern at 1 30, and the technician does not reinitialize automode at 1 34, the 
ROI is determined by the technician at 129, and the processing above is repeated. If a match was 
not foxmd at 126, and more than a small number of frames have failed to result in a matching 
pattern at 130, and the technician does reinitialize automode at 134, the automode is returned to. 
If a match was not found at 126, and more than a small nvunber of frames have not failed to result 
in a matching pattem at 130, the ROI is determined to be the ROI from the previous fame at 132, 
the technician may reinitialize automode at 1 36 and return to automode, or the technician may not 
reinitialize automode at 136, and a frame is selected for analysis at 122. And the processing 
described above is repeated. 

INDUSTRIAL APPLICABILITY 
[0061] The present invention finds industrial applicability in the field of automated 
determination of a Region Of Interest (ROI) within an Image. More particularly, the present 
invention finds industrial applicability in the field of automated determination of an ROI of an 
ultra sound image of a carotid artery. 

SCOPE OF THE INVENTION 
[0062] While the invention herein disclosed has been described by means of specific 
embodiments and applications thereof, numerous modifications and variations could be made 
thereto by those skilled in the art without departing from the scope of the invention set forth in the 
claims. 
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